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ABSTRACT 



■ Under the assumption that accretion onto massive black holes powers active galac- 

Cn| \ tic nuclei ( AGN) , the mass function (MF) of the BHs responsible for their past activity 

is estimated. For this, we take into account not only the activity related to the opti- 
cally selected AGN, but also that required to produce the Hard X-Ray Background 
(HRXB). The MF of the Massive Dark Objects (MDOs) in nearby quiescent galaxies 
| is computed by means of the most recent results on their demography. The two mass 

functions match well under the assumption that the activity is concentrated in a single 
significant burst with A = L/ LEdd being a weakly increasing function of luminosity. 
This behaviour may be indicative of some level of recurrence and/or of accretion rates 
insufficient to maintain the Eddington rates in low luminosity/low redshift objects. 
Our results support the scenario in which the early phase of intense nuclear activity 
occurred mainly in early type galaxies (E/SO) during the relatively short period in 
which they had still an abundant interstellar medium. Only recently, with the de- 
cline of the QSO luminosities, did the activity in late type galaxies (Sa/Sab) become 
statistically significant. 



1 INTRODUCTION 



S-h ' The number counts of AGN's and the intensity of the backgrounds at high energies show that the activity in nuclei of galaxies 
5^ i was much higher in the past than in the local universe. If we accept the paradigm that nuclear activity in galaxies is sustained 
by accretion onto massive BHs (see Rees 1996 for a review), then the problem of the location and discovery of the remnants 
of such past activity immediately arises. There is evidence for the presence of Massive Dark Objects (MDOs) in the centers 
of most, if not all, nearby galaxies with a large spheroidal component. The mass of MDOs has been evaluated by using very 
high resolution spectroscopy and photometry of the centers of nearby host galaxies. Observations with HST have allowed a 
significant breakthrough in the angular resolution and have led to an enormous increase in sensitivity and precision of mass 
estimates (Magorrian et al. 1998, hereafter M98; van der Marel 1997, hereafter VdM97; Ford et al 1997, hereafter F97; for 
a review, see Kormendy and Richstone 1995, hereafter KR95). As a result, a large number of MDOs have been detected 
suggesting that we are discovering the fossil of the past nuclear activity. Here we will assume that MDOs are the BHs now 
dormant after a shining past. 

The main purpose of this paper is to show that the mass function of the BHs we infer from the observations of the past 
activity of AGN and QSOs corresponds to the MDO mass function of the local non active galaxies. 

M98 have successfully exploited the very high resolution of HST photometry and ground based spectroscopy of 36 E and 
SO galaxies in order to estimate their MDO masses. In spite of the large scatter in the data, they also confirmed the existence 
of correlation between the mass of the hot galactic component M 3p h and the MDO/BH mass, Mm do, already suggested by 
Kormendy (1993) and KR95 and successively claimed by M98 and van der Marel (1998) (hereafter VdM98, by analysing a 
sample of 46 early type galaxies). These large samples, together with several smaller ones (e.g. F97; VdM97; H98) allow to 
estimate the distribution function of the ratio Mm do /M Bp h and then to evaluate the mass function of the MDOs (see section 
2.1). 
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A different approach to evaluate the MDO/BH mass function relies on the hypothesis that radio emission from the nuclei 
of radio quiet galaxies is related to the mass of their MDOs. As a matter of fact, a strong correlation between MDO masses 
and nuclear radio luminosities has been found by Franceschini et al (1998). This correlation, in connection with the radio LF 
of the nuclear emission of radio quiet galaxies, can be used to probe the MDO mass function (see section 2.2). 

In section 3 we will determine the mass function of the material accreted onto BHs by exploiting the knowledge on the 
evolution of AGN/QSO LFs. Reliable luminosity functions and cosmic evolutions are presently available for optically and 
soft X-ray selected objects. On the other hand, we must also consider a class of heavily absorbed AGN, under-represented in 
optical and soft X-ray surveys, that shows up in hard X-ray surveys. These absorbed AGN are the most likely contributors of a 
major portion of the intensity of the 2-50 keV X-ray background (HXRB) through BH accretion energy output. In fact, optical 
identifications of serendipitous sources detected by BeppoSAX in the 5-14 keV band show that absorbed AGN significantly 
contribute to the HXRB (Fiore et al. 1998). Moreover, independent observational evidence supports unified schemes for the 
AGN, whereby a large fraction of absorbed (type II) AGN is expected (see Antonucci 1993 for a comprehensive review). Thus, 
we will derive the mass function of relic BHs by including not only the mass accreted on type I or optically selected AGN's, 
but also that accreted on type II AGN's. 

In section 4 the comparison of the MF of the dormant BHs to the accreted MF will be used to cast light on the 
characteristics of the evolution of the nuclear activity in galaxies of different morphological type. 

We will adopt H = 70km/s Mpc" 1 unless otherwise stated. Moreover, ft = H / in units of 100 kms~ 1 Mpc~ 1 . 



2 ESTIMATES OF THE MDO MASS FUNCTION 

2.1 From the optical luminosity function to the MDO mass function 

Observations and subsequent analysis have shown that MDOs are quite common in galaxies with significant spheroidal 
components and that their masses are correlated with the spheroid masses, though with large scatter (Kormendy 1993, KR95; 
VdM97; VdM98; M98; F98). The mean value of M M Do/M sph is still under debate and it ranges from (M M Do/M sph ) ~ 10~ 2 
(M98) to (MMDo/M S ph) ~ 2 x 10~ 3 (H98). The uncertainty is related to different assumptions on dynamical models, 
particularly on the two-integral phase-space distribution function (VdM97; H98; M98). Exploiting very high resolution HST 
photometry and ground based spectroscopy, M98 estimated the mass of 36 MDOs in the centers of nearby galaxies, mainly 
E and SOs. They found that the distribution of the ratio x — Mmdo IM sph can be described by a Gaussian distribution in 
log x: 



where TV is a normalization constant. The best fit values found by M98 are log x = log{ Mmdo /M ap h) = —2.28 and a = 0.5. 
Analyzing the very high resolution photometry of 46 early type galaxies and reproducing it by means of the models proposed 
by Young (1980), VdM98 found log Mmdo — —1.83 + log Ly with rms scatter a ~ 0.3 dex. Using standard values of 
M/Lv — 5.5 h(Lv / L*) 0,25 for spheroids (see below), this corresponds to (log MMDo/M sp h) — —2.64. 

We have constructed a sample of 30 galaxies, for which at least two independent mass estimates of their MDOs are 
available in the literature (KR95, VdM97, VdM98, H98, F98, M98). Taking the geometric mean as the best estimate of the 
ratio x, we have found that a Gaussian distribution in log x, with {log x) = —2.60 and o — 0.3 is a good description of 
the data. This result, quite close to that obtained by VdM98 and in agreement with the estimates of H98 and F98, will be 
adopted as a reference distribution. 

The estimate of the spheroid mass function in the local universe is a preliminary step in the evaluation of the MDO MF. 
We determine the luminosity function of the spheroids by using the total- luminosity function (L to t) and the estimates of the 
ratio of the spheroidal to total luminosity L sp h/L to t for different types of galaxies (E, SO, Sa/Sab, Sbc/Scd). The LFs of early 
type and spiral galaxies are well known in the range —23 < Mb < —17. In fact, results of different groups (e.g. Efstathiou 
et al., 1988; Loveday et al., 1992; Zucca et al., 1997; Geller et al., 1997, Heyl et al. 1997), based on different galaxy samples, 
well agree in the luminosity range relevant to our purpose. In order to represent the LFs of the various Hubble types, we use 
standard Schecther functions whose parameters L* and a are given in Table 1 for ft = 1 
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Figure 1. Left panel: a) The mass function of the spheroids. The solid line is the total MF, the short-dotted, dot-dashed, dashed and 
dot-long dashed lines refer to spheroids in E/SO/Sab/Sbd. Right panel: b) The MDO mass function, (froMFdlog Mjj: i) derived from 
eq (1) with the values of (x) and <r given in section 2.1 (solid line), ii) relative to the contribution of MDO's in E/SO's (long dashed 
line) Hi) derived by Franccschini ct al (1998) (solid circles) iv) derived by M98 (dotted line). 

The ratio L sp h/L t ot has been investigated by several authors (Simien and de Vaucouleurs 1986, Kent 1985, Kodaira et 
al. 1986) and it has been found to vary along the Hubble sequence. There is a reasonable agreement on the average values, 
though intrinsic scatter is significant, particularly in the late types. With the values reported in Table 1, the total luminosity 
density in spheroids turns out to be ~ 30% of the total luminosity density, in agreement with Schechter and Dressier (1987) 
and Fukugita et al. (1998). 

From the spheroid LF we infer the MF by adopting M sph /Lv ^ 5.5 h(Lv / L+) 0,25 , corresponding to M sph /L_g ~ 
6.9 fe(L s /L*) - 25 with (B-V)=0.90. This law is compatible at 2 a level with (M sph /L v ) oc L^ 18±0 ' 03 found by M98 and 
(M sph /L B ) oc L° B 35±0 - 05 claimed by van der Marel (1991); consistent with stellar evolution models. The spheroid MF is 
presented in figure la. The overall shape of the MF is relatively flat at M sp h ^ 10 11 Mq and exhibits an exponential decline 

for M sph > 10 n M Q . At M sph > 5 x 10 10 M© the mass function is dominated by spheroids in E and SO galaxies. The Sa/Sab 
galaxies exhibit an exponential decline of their MF for M sph > 10 10 M Q and Sbc/Scd galaxies for M sph > 10 8 M Q . 

The MDOs mass function (hereafter OMF) has been computed by convolving the MF of the spheroids with the adopted 
log Gaussian distribution with (Mmdo /M sp h) = —2.6 and a — 0.3. As it is apparent in figure (lb), the scatter of the ratios 
x = Mmdo /M ap h softens the exponential fall off of the luminosity function; the broader the distribution the gentler is the 
decline. Of course, the convolution relies on the assumption that scatter reflects a real complexity of the physical processes 
leading to the formation of the BHs in galaxy centres. 

The total mass density amounts to pmdo — 8.2 x 10 5 Mq/Mpc 3 h 2 , with large fraction ~ 75% due to MDOs in E and 
SO galaxies. Of the remaining 20% is due to MDOs in Sa/Sab galaxies and 5% to MDOs in late type spirals. The results 
for MDOs in spirals are consistent with the upper limits that Salucci et al. (1998) have derived from the analysis of a large 
number of high quality rotation curves. 

The differences with the MFs derived using the distributions proposed by M98 are apparent from figure lb. The M98 
log Gaussian law predicts a large number density of BHs with mass larger than 10 10 Mq and yields a total mass density 
Pmdo ~ 3 x 10 6 Mq/Mpc 3 h 2 , which is higher than the estimate of the local mass density of the accreted matter, based on 
counts of optically selected AGN, by a factor of about 10-15 (see section 3). 

A contrasting result has been obtained by Franceschini et al. (1998), whose mass function, estimated from the optical 
LF of E/SO galaxies and the correlation between Mmdo and the luminosity of the bulge of the host galaxies, is substantially 
below our determination and yields a lower total mass density pmdo ~ 10 5 Mq/Mpc 3 h 2 . This is partly due to the adopted 
LF of E and SO galaxies and partly to the fact that they neglected the broadness of the distribution of the ratio Mmdo / L sp h- 

2.2 From the radio luminosity function to the MDO/BH mass function 

A different but complementary approach to the evaluation of the MF of BHs is based on the idea that the radio emission from 
cores of normal galaxies is related to the mass of the hosted inactive BHs. Actually, accretion onto dormant BHs occuring at 
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very low rates is possibly advection-dominated with low radiative efficiency (Rees et al 1982; Abramowicz et al 1988; Narayan 
and Yi 1995 a,b). In this regime, a hot ion torus (T ~ 10 12 ° K) forms, while the electrons, coupled to the ions only through 
two-body Coulomb effects, remain at much lower temperature (T ~ 10 9 ° K) and emit at radio frequencies via thermal 
synchrotron. Fabian and Rees (1995) pointed out that the thermal synchrotron emission of the electrons might be detectable 
in cores of quiescent galaxies. They also noted that most of the elliptical and SO galaxies exhibit a low-power emitting radio 
core. The ADAF model predicts for the core radio power P v oc v 2 ^ 5 M 4 ^ 5 (e.g. Mahadevan 1997). The dependence on 
frequency is quite close to the relationship P„ oc v 1 ^ 3 found by Slee et al (1994) for a sample of radio cores of E/SO galaxies. 
However the ADAF model should be taken with some caution, since new high resolution radio and sub-millimeter observations 
of three giant elliptical galaxies significantly disagree with its canonical predictions, although the possibility of explaining the 
observed spectra with modifications of the canonical ADAF is not ruled out (Di Matteo et al. 1998). With Bondi accretion 
rate M oc M% H p(oo) / 'c 3 a {oo) the ADAFs yield radio powers P oc Mg%, if the density and sound velocity at the boundaries 
of the accretion flow are independent of the mass of the spheroid. When the latter depend on spheroid mass, p oc M sp h 
and c s oc oc M^^-, we have P oc Mgjj. Given that self absorbed processes (e.g. jets) yield P oc M% H , a relationship 

P oc Mg H with 2.0 ^ ft ^ 2.2 is expected under rather general conditions. 

Franceschini et al. (1998) find a significant correlation between the radio power of the cores P CO re and the estimated BH 
mass in 8 objects with P CO re oc MgJ^ ' 6 . In a new analysis, we searched for high angular resolution radio observations of the 
galaxies with at least two MDO mass estimates, as specified in the section 2.1. The radio observations reported by Wrobel 
(1991), Wrobel & Heeschen (1991), Slee et al. (1994), Sadler et al. (1995) allow the selection of a sample of 15 objects with 
detections both in mass and in radio emission and 6 objects with upper limits in the latter. 

The correlation and the linear regressions have been computed with the package ASURV (LaValley, Isobe & Feigelson 
1992). The probability of correlation by chance turns out to be P ^ 0.4%. The best estimate of the regression is: 

log PlZT - 19.5 + 2.2l :£ log (^-) , (2) 

with a slope close to the theoretical inferences for the exponent a (figure) 2a. It is worth noticing that the same value of the 
slope has been claimed also by McLure et al. (1998). 

In order to obtain the MF of BHs in nearby galaxies, we need to convolve equation (2) with the radio luminosity function 
of the nuclei of quiescent spheroidal galaxies. Sadler et al. (1988) surveyed 114 nearby E and SO radio-quiet galaxies at 5 GHz 
to search for weak radio sources and computed the total power function 4>{Pt)dlog Pt at 5 GHz down to log Pt(W/Hz) = 19.6 . 
More recently Slee et al. (1994) observed a large fraction of objects from the same sample at an higher angular resolution; 
this has provided the Pt vs. P CO re relationship: log P CO re = 19.5 + 0.78 (log P t — 19.5) for log Pt > 19.5 and log P CO re = log Pt 
at lower powers. We can then convert the radio total power function of Sadler et al. (1988) into the core power funtion (PF) 
(figure 2b) suitable for the present study. 

By using equation (2), we pass from the radio core PF to the MDO/BH mass function (hereafter RMF). The result is 
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Figure 3. Left Panel: a) The elliptical's OMF mass function compared to the mass function, ipRMF dlog Mbh, derived from the E/S0 
cores radio PF (filled circles). Also shown the RLF derived by Franceschini et al. (1998) Right Panel: b) The total mass function of the 
relic BHs, <f>AMF dlog Mjg, derived from the past activity of AGN's (thick dashed line); also shown the contribution from the AGN's 
originating the HXRB (thin-dashed line) and the total mass function for h = 1/2 (dotted line). 



shown in figure (3a) alongside with its uncertainties. Let us notice that for obtaining the BH mass function, the available 
statistics can only roughly estimate the uncertainties related to systematic errors in the procedure. While the errors in density 
represent only the statistical errors of the radio LF, we estimate the errors in mass from the relationship log Mbh versus 
log Pt, taking into account its uncertainties in the normalization and in the slope. The errors turn out to be about ±0.2 dex. It 
is worth noticing that the core radio luminosity function exhibits a gentle slope and no exponential decline at least for powers 
P ^ 10 24 W/Hz. This ensures that the possible scatter in the log P cor e versus log Mbh relationship does not significantly 
affect our results. 

Franceschini et al. (1998) estimated the MF using the total radio luminosity function of the E and SO galaxies and the 
correlation of the total power with the BH mass Pt oc M B 'f^- The resulting MF stays significantly below our estimate for 
M B h > 2 x 10 8 M Q and the predicted total mass density is smaller by a factor of about 3, mainly because of the steeper 
adopted relationship P t — M B h ■ 

In figure 3a we also plot the MF determined via the optical LF. The comparison should be done with the lower curve, 
which refers to the E/S0 galaxies as the RMF. Although the two MFs are not completely independent, nonetheless, their 
agreement supports the reliability of both estimates. The dependence of the two methods on the distance scale is quite similar 
and thus equivalent matches can be found also for h = 0.5. It is worth noticing that the OMFs derived by using the distribution 
of log x proposed by M98 (see eq. 1-2 and figure lb) cannot be reconciled with the RMF. 



3 THE LOCAL MASS FUNCTION OF THE DORMANT BLACK HOLES 

The mass density associated to luminosities higher than L is : 

k f 1 + 2 f'Lmax 

p(> L) = Jk J (i+,)3(i + » 2) i/ 2 d2 J L Ln ^ z ^ dL ( 3 ) 

where L is the luminosity in a certain band and kt> i is the corresponding bolometric correction, O is the density parameter 
and n(L, z) is the luminosity function at redshift z. As for the conversion efficiency of rest mass of the accreted matter into 
energy, we adopt e = 0.1, unless otherwise stated. Putting L = L m i n ~ 10 44 erg/s (e.g Pei, 1995), we get the local total 
mass density pbh- As shown by Soltan (1982) this quantity can be written in term of source counts and is independent of 
the cosmological parameters. 

The mass function of the matter accreted onto AGN during their activity can now be derived by means of the large 
amount of data available on spectra and on the evolution of the luminosity function. The increasing evidence of the presence 
of MDOs (i.e. inactive BHs) in a large number of galaxies with significant spheroids favors the scenario whereby the nuclear 
activity is a relatively short phase occurring in a large fraction of galaxies. The continuity paradigm is then excluded. 
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In the hypothesis of single short events it is likely that the AGN are observed in the highest state of activity, then, we 
can evaluate the mass of the BH associated to a luminosity L, (hereafter the bolometric luminosity) through the assumption 
L — \L E dd and L Edd = M B hc 2 /t B - The constraint for the latter can be obtained as it follows: if we assume A constant during 
the accretion, we get: 

M BH (t obs ) = M BH (t m )o X p (^g), (4) 
and 

L{t ob3 ) = ^—M B H{t bs) = ^—M BH (t ln )exp (5) 
Ze tE \et E / 

where t B = (arc) / (4-Km p ) ~ 4 x 10 8 yr, and t in and t ob3 are the initial time of the bright phase and the observing time 
respectively, and r = t b s — tin- Then, by requiring M B H{t bs) S> M B H{Un) we get the constraint: r > (2etfi)/A. 
Following equation (4), we write the QSO/AGN luminosity function ip(L') as 

' L ' J„IJ^ T\ r l °3 M'a(L') 

(6) 



/ VK U L= I M a (j>AMF{M a )dlog M a 

J J — oc 



where 4>amf{Mo) is the mass function of the accreted mass, M a = M B h- Since A is likely a weakly increasing function 
of the luminosity, let us then set: 

\(L)=L/L E = W jil ° 9 L ~ 49) ; 7^0.2 (7) 

that follows the findings that the most luminous QSOs radiate at about the Eddington limit, while low luminosity AGN, 
L ~ 10 44 erg/s, show A = L/L E ~ 0.1 - 0.05 (see Padovani (1989) and Wandel (1998)). Then, from equation (6) and (7) the 
MF of the relic BHs can be written as: 
j.,*, nji Zn(10) ip(> L) dlog (ip(> L)) „ 

mBH )dlog M BH = {1 { _^ K M J g dl ( J g \ L) » dlog M BH (8) 

Next step is to insert into equations (3) and (8) the information on the luminosity function and its evolution. For the 
optically selected AGN we adopt the luminosity function and its cosmic evolution for the B-band given by Pei (1995), which 
is well defined at least up to z ~ 3.5. The bolometric correction for the B-band has been taken kboi{B) = 13, on the basis of 
the spectra reported by Elvis et al (1994). The total mass density turns out to be p BH — 2.0 x 10 5 M Q /Mpc 3 , close to the 
estimate of Chokshi & Turner (1992), who however adopted k bo i(B) — 16.5. 

The mass density of matter accreted on massive BHs powering optically selected AGN has been often used as a reference 
for the total mass density in BHs. However, as pointed out by Granato et al. (1997), the total amount of the matter accreted 
onto AGN should include that required to produce the 2-50 keV background (HXRB). In fact, the surveys in the optical and 
in the soft X-ray bands tend to select expecially type I AGN and QSOs (Hasinger 1998) and to lose a significant fraction of 
absorbed active objects, the so called type II AGN. 

Setti and Voltjer (1989) suggested that these AGN yield most of the HXRB (let us remind that type I AGN are well 
known to be minor contributors to the HXRB). A direct evidence comes from Fiore et al (1998), that, in a survey with 
BeppoSAX in 5-14 keV band, found about 150 sources corresponding to a surface density of 20 objects/sq degree at the flux 
limit F ~ 6 x 10~ 14 erg/s/cm 2 . So far, there have been optically identified only 9 sources, 5 of which turn out to be QSO's 
and 4 "type II" AGN. The surface density implied by these detections can explain about 40% of the HXRB. 

A suitable estimate of the mass density underlying the intensity of the HXRB between 2-50 keV can be obtained through 
the relationship 

Ebol _ k bo l (1 + Z e ) I 

k(z e 

The intensity in the 2-50 keV band is I ~ 1.9 x 10~ 6 erg cm~ 2 sec -1 , k bo i ~ 14.5 is the bolometric correction for 2-50 keV 
band (see Elvis et al. 1994; Bassani et al. 1998), z e is the effective emission redshift, which can be assumed z e ~ 1 and 
k(z e ) is the k-correction for z — z e . The total associated mass density amounts to p B H ~ 3 — 5 x 10 5 Mq/M^c -3 , where the 
uncertainty reflects the uncertainty on z e , which disappears once n(L, z) is known. 

Direct information on n(L, z) for objects selected in hard X-ray bands is still scarce. However, assuming that the HXRB 
is chiefly ascribed to type II AGN and that the unified schemes are basically valid, then the soft X-ray observations and the 
shape of the HRXB significantly constrain the luminosity function and the evolution of the type II AGN. In order to compute 
the mass function, we use the models proposed by Comastri et al. (1995) and by Celotti et al (1995) to reproduce the HXRB. 
Comastri et al (1995) used a LF in the 0.3-3.5 keV band, while Celotti et al (1995) one referring to the 2-10 keV band. We 
adopted kbol = 25 and kbol = 38 respectively, on the basis of the spectra presented by Elvis et al. (1994). Both models assume 
luminosity evolution L(z) cx (1 + z) k with k ~ 2.6 — 3.0 and a similar redshift cutoff z ~ 2 — 3 beyond which the evolution 
stops. The constraints imposed by the shape of the HXRB and by the LF and the evolution of the soft X-ray selected AGN, 



J^bot ri-bol y± T 1 /n s 

PBH ~ ~ — ~ j j / -■ 9 
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Figure 4. Left panel: a) The AMF derived by assuming A = 1 (dotted line) and A = 0.2 (dashed line). In both panels the solid line and 
filled circles are defined as in fig. (3a). Right panel: 4b) The AMF derived according to eq. (7) compared with the MF of the relic BHs 
by Cavaliere and Vittorini (1998) (dotted line). 

compel the redshift distribution to be quite similar for both models. Let us notice that the shape of the mass function and 
the mass density pbh(X — ray) ~ 4.2 x 10 5 M Q /Mpc 3 results comfortably very similar in both models. 

In figure 3b we show the mass function derived for A = 1 and A given by eq (7). The corresponding total mass density is 
Pbh — 6.5 x 10 5 Mo/Mpc 3 . It is apparent that optically selected AGN dominate at Mbh > 10 9 Mg, while BHs associated to 
the X-ray absorbed AGN dominate the MF at M < 10 8 M Q . Unlike the total mass density, the mass function of the accreted 
matter (hereafter AMF) does depend on the cosmological parameters. However, as shown it figure (3b), if we decrease h down 
to 0.5, the AMF, for Mbh > 5 x 10 9 , changes slightly . 

The local LF of AGN in the X-ray bands do not extend below sim4 x 10 42 erg/s so that our estimate includes only 
objects with L > 10 44 erg/s and with Mbh > W 7 M Q , assuming A ~ 0.1 at low luminosities. In order to probe masses below 
this limit we should include the, presently not very well-known contribution from low luminosity active nuclei, such as the 
LINERS. 

In figure (4a) we compare the MF derived from the past activity, assuming A = 1 and A = 0.2, to the mass functions 
derived from the optical and radio LFs. The slopes are quite different and the agreement may be considered acceptable only 
in the range 10 8 Mq — 10 10 Mq. This suggests that the assumption of a constant A = L/Le is not adequate. 

Instead, by using eq. (7), which is motivated on observational grounds, the agreement is excellent. In detail, since the 
total mass density is constant, the individual masses of low luminosity objects increase while their number decreases, and the 
AMF of the inactive BHs turns out in good agreement with the MF of the local MDOs (see figure (4b)). This dependence 
may be suggestive of recurrent activity confined to low mass objects. 

4 DISCUSSION AND CONCLUSIONS 

The agreement found between the mass functions derived from investigations on MDOs resident in local galaxies and the 
mass function of the BHs inferred from the past activity of AGN is based on very simple and sound hypotheses: i) the nuclear 
activity is a single short event; ii) the spectra of the AGN do not greatly depend on redshift; Hi) the mass-radiation conversion 
efficiency of accretion e ~ 0.1; iv) the HXRB is produced by absorbed AGN; v) A = L/Le is a weak increasing function of 
the luminosity. 

The inclusion in the estimate of the mass deposited in BHs by the activity related to the HXRB is mandatory, indepen- 
dently of the specific model of HXRB one adopts. In fact, the optically selected objects contribute a minor fraction £ 20% 

of the HXRB, whose observed intensity implies pbh (X - ray) - 3 - 5 x 10 5 Mo/Mpc" 3 . 

Let us comment that the idea that ADAF accretion would be responsible for the HXRB (Di Matteo and Fabian, 1997, 
see also Haenhelt et al 1998), which implies a quite larger pbh(X — ray) and a different AMF with respect to the one we have 
derived, runs against the evidence from the new hard X-ray surveys that have found that obscured objects are responsible 
for the emission of at least 20 % of the HXBR, yet at a rather bright flux limit (Fiore 1998). 
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The MF computed using the distribution of the ratio Mbh /M sph obtained by M98 is in contrast with that derived from 
the correlation found between MDO mass and radio core power. If we neglect this fact, the OMF predicted by means of the 
log normal distribution of the Mbh /M sp h ratio of M98 can match the AMF by assuming an efficiency e ~ 0.02, a factor 5 
less than the standard value. Alternatively, the bolometric corrections should be increased by the same unlikely factor. 

Large mass densities could be explained by assuming that the mass accreted during bright phases that shows up directly 
in the optical counts or in the integrated HXRB is only a small fraction of the total mass of the BHs; in this case most of 
the mass must be accreted in a silent phase (in the optical and in the X-ray bands) .However, If the silent phase occurred 
before the bright one, the problem of BH formation at early time would become more intricate and A << 1 is required at all 
luminosities. If it is postponed, then extreme obscuration is required. 

The dependence on H a of the mass functions is interesting. The shape of the AMF depends slightly on H a while the 
corresponding total mass density of the accreted matter is independent of H a . On the other hand, the mass functions derived 
from the local radio power function and optical LFs strongly depend on H a , with Mmdo ^ 1 and $ <x /i 3 while the mass 
density in dormant BHs pBH oc h 2 . With the present uncertainties in the derivation of the mass functions it is not possible 
to discriminate among different values of the Hubble constant. In order to counterbalance the decrease by a factor of about 
2 of the total luminosity density due to the change from H a = 70krn/s Mpc~ x to H a = 50km/ s Mpc" 1 , the average ratio 
Mbh /M sp h must be increased by the same factor. Thus the match between OMF and AMF is obtained with log x a = —2.35 
and a = 0.3 in eq (1). The agreement with the RMF is obtained assuming log Plo" z - 19.2 + 2.0 log (M B h/W s M ), which 
is a good fit to the data. 

However it is conceivable that in the near future additional high resolution photometric, spectroscopic and radio observa- 
tions and very deep hard X-ray surveys will allow to significantly reduce the uncertainties in the estimate of OMF, RMF and 
AMF. At that point the possibility of balancing changes of the distance scale acting on some parameters will be interestingly 
restricted. 

The activity patterns have been investigated by several authors (e.g. Cavaliere and Padovani 1988; Small and Blandford 
1992; Haehnelt and Rees 1993; Cavaliere and Vittorini 1998). Cavaliere and Vittorini (1998) propose a scenario in which 
the energy output, at early times are governed by hierarchically growing environment, while later, the fall in the average 
luminosity is related to intermittent accretion governed by galaxy-galaxy interactions. Their prediction of the local MF, 
reported in figure 4b, is below our estimate in the mass range 10 7 Mq £ Mbh ^ W 9 Mq. The derived total mass density is 

Pbh — 3.2 x 10 5 Mq Mpc~ 3 , falls short by a factor two the total mass accreted onto type I and type II AGN. The difference 
may be due to the fact that BHs in early type spirals are not considered by Cavaliere & Vittorini (1998). 

In conclusion our analysis strongly supports that the nuclear activity is in general a single event, with bolometric lu- 
minosities close to the Eddington luminosities A = L/L E dd but increasing increase from A ~ 0.1 for L ~ 10 44 erg sec -1 to 
A ~ 1 for L > 10 48 erg sec -1 , as found also by independent studies (see e.g. Padovani 1989). This suggests that recurrent 

activity occurs in low luminosity objects and/or that their accretion flows are supply-limited and cannot sustain Eddington 
luminosities. 

In our framework type II AGN of small and moderate mass are responsible for most of the HRXB. These objects can 
reside in early type galaxies but also in spirals with still significant spheroidal component (Sa/Sab). On the other hand, bright 
objects such as QSOs at significant redshift should be hosted preferentially by E/S0 galaxies (or by their precursors). This 
has been recently confirmed by a study of the host galaxies of QSOs with Mr > —24 (McLure et al. 1998). These author 
found that a large fraction of the host galaxies are massive elliptical galaxies. This facts add complexity to the problem of 
the AGN evolution. Not only coordination is required to mimic the luminosity evolution, but also a kind of coordination is 
required to pass from the high luminosity, high redshift optically selected AGN, preferentially hosted in E/S0 galaxies, to the 
low luminosity objects , such as the local Seyfert galaxies, hosted preferentially in Sa/Sab galaxies. The discussion of these 
aspects in relation with the problem of BH formation will be presented in a subsequent paper (Monaco et al 1998). 
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